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Pure glue QCD is formulated on a 2+1 dimensional trans- 
verse lattice, using discrete light-front quantization. The 
transverse component of the gauge fields is taken to be com- 
pact, but in a linearized approximation with an effective po- 
tential. The (rest-frame) QQ potential is evaluated numeri- 
cally using Lanzcos matrix diagonalization. We first discuss 
the strong coupling limit analytically and then present nu- 
merical results beyond the strong coupling limit. The phys- 
ical origin of confinement on a transverse lattice depends on 
the orientation of the external charges: For longitudinally 
separated charges, confinement arises from the instantaneous 
Coulomb interaction and for transversely separated quarks a 
string of link-fields forms. In the general case one obtains 
a superposition of both effects. Despite the asymetry in the 
microscopic mechanism, already a very simple ansatz for the 
effective link-field potential provides an almost rotationally 
invariant QQ potential. The momentum carried by the glue 
depends strongly on the orientation of the external charges, 
which might have observable consequences. 



I. INTRODUCTION 

Light-cone coordinates |l| are the natural coordinates 
for describing high-energy scattering . The immense 
wealth of data on nucleon and nucleus structure functions 
thus strongly motivates to understand QCD on the light- 
cone. The transverse lattice formulation of QCD |^-||] 
is a particularly promising approach towards this goal. 
Among the most appealing features of this approach to 
light-front QCD is that confinement emerges naturally in 
the limit of large lattice spacing HJl^ . 
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FIG. 1. Space time view of a transverse lattice 

In the transverse lattice formulation, the transverse 
space directions are discretized, while the longitudinal 
[i.e. the = (a;° ±x^)/\/2] directions are kept continu- 



ous (Fig.|l|). While this seems to be a natural procedure 
when quantizing on the light-front, the procedure is obvi- 
ously not manifestly rotationally invariant and one might 
ask oneself whether rotationally invariance is recovered in 
the continuum limit This issue becomes even more 
relevant, when one considers "practical" calculations (in 
contrast to the infinitely complicated continuum limit), 
i.e. calculations where the lattice spacing is not necessar- 
ily infinitesimally small because of numerical limitations 
or where one makes approximations on top of the dis- 
cretization. 

In this paper, we will consider one specific observable, 
namely the rest-frame potential energy of an (infinitely 
heavy) QQ pair coupled to the gluon system. As has 
been shown in Ref . , this observable can be extracted 
from a light-front Hamiltonian by considering a QQ pair 
which moves with uniform velocity and where the sepa- 
ration between the quark and the antiquark is kept fixed 
(the corresponding formalism is briefiy summarized in 
the Appendix). 

After setting the general formalism, we will proceed 
to calculate the QQ potential away from the strong cou- 
pling limit. The main issues there are to find out whether 
linear confinement persists and whether rotationally in- 
variance for the potential is being restored. 



II. THE HAMILTONIAN 

The light-front Hamiltonian for compact QCD on a 
transverse lattice has been introduced in Ref. For 
pure glue QCD in 2-1-1 dimensions one finds 

= / dx~ / dy~ : tr [J„(x") \x~ - y~\ J„(y")] : 

+ Veff{U), (2.1) 

where 



(2.2) 



and Un are the link fields, which are quantized matrix 
fields and satisfy the usual commutation relations. Ide- 
ally, one would like to work with [/„ G SU{N), but in 
practice this is very complicated |^ so in practice one 
may prefer to work with an unconstrained complex ma- 
trix field and instead add an effective constraint term 
Veff{U) to the light-front Hamiltonian. In the case of 



1 



N — > oo, in the classical limit, Veff{U) can be taken of 
the form 

V,'^f{U) = C2^tr : [UlU„] : +C4 ^ tr : [C/tC/„[/tC/„] :, 

n n 

(2.3) 

where C2 = — 2c4 and C4 — > c», which provides an effec- 
tive potential that is minimized for [/„ e U{N). In the 
N ^ 00 limit, the dif fere nce between U{N) and SU{N) 
is irrelevant and Eq.( ^.3[ ) is thus suitable for enforcing 
the SU (N) constraint in the classical limit. One might be 
tempted to try a similar ansatz for the LF quantized case. 
There are several reasons why a different form for the ef- 
fective potential is more useful. First, if one still attempts 
to work with an effective potential of the above form then 
physical states would necessarily look extremely compli- 
cated, which has to do with the fact that the above ansatz 
for the effective potential corresponds to working clo se to 
the continuum limit. Thus even if the ansatz in Eq.( |2.3| ) 
would work in principle, it would most likely not be very 
practical. 

However, since a "Mexican hat" potential corresponds 
to a situation where one is working with the false vacuum, 
it is questionable whether a physical situation where a 
particle runs at the bottom of a Mexican hat can be 
described at all by a LF Hamiltonian, using degrees of 
freedom expanded around the origin. 

For these two reasons, it makes more sense not to con- 
sider Un as the bare link field, but instead think of it as 
some kind of blocked or smeared variable. The blocking 
has several consequences. First, the SU{N) constraint 
gets relaxed which reflects itself in the fact that the ef- 
fective potential is no longer just a narrow valley p2[ . 
Secondly, using smeared variables, it might be easier to 
cover large physical distances with only few degrees of 
freedom. The price one has to pay for these advan- 
tages is that the effective potential gets more complex 
and in general more terms are necessary than shown in 
Eq.(2.3). In Ref. an attempt has been made to fit the 
effective potential to the glueball spectrum by making an 
ansatz which includes all operators up to dimension four. 
Since this work is a first study of the rest-frame QQ po- 
tential and, as we will discuss below, the QQ potential 
turns out to be rather insensitive to terms of dimension 
greater than two in Ve//(J7), we will instead only con- 
sider a much simpler ansatz in the following and keep 
only the quadratic term 



the gauge field /7„ will be expanded in of a plane wave 
basis of complex matrix fields. For simplicity, zero mode 
degrees of freedom will be omitted in this procedure, but 
we impose a color singlet constraint at each site 



Veff{U)^C2tr{U^U) 



(2.4) 



III. NUMERICAL PROCEDURE AND RESULTS 

Once one has made an ansatz for the Hamiltonian, one 
needs to find approximate solutions. In this work, we 
will employ DLCQ H for this purpose, i.e. on each link 



Qn\phys) = dx Jn{x )\phys) = 0, 



(3.1) 



which is known to results after eliminating zero mode 
degrees of freedom in the limit of a large longitudinal 
interval [| In the large N limit, the color singlet 

constraint is easily satisfied by using a basis of states 
that can be written as traces over color indices. This has 
been shown in Ref. |^ and will be used in this paper. 

To some extend, one can justify the omission of zero 
modes in the sense of Refs. p"6[-p]| by formulating the 
theory in terms of an effective light-front Hamiltonian, 
where high-energy degrees of freedom have been inte- 
grated out. While this procedure certainly takes care of 
some of the zero-mode dynamics, it is not clear whether 
this allows one to completely omit zero-modes as dynam- 
ical degrees of freedom. Therefore, one should regard the 
omission of zero-modes (as dynamical degrees of freedom) 
in this paper as an ad hoc approximation which is applied 
to make the problem simpler. Including dynamical zero- 
modes is one of the many possible improvements that one 
could consider as an extension of this work |22|. 

In principle, there can be an arbitrary number of gluon 
quanta on each of the links of the transverse lattice. The 
only constraint is abovementioned color singlet require- 
ment on each site. In order to simplify the numerical 
calculation, we will first calculate the QQ potential un- 
der the (ad hoc) approximation that each link contains at 
most one gluon quantum. In combination with the color 
singlet requirement this implies that a QQ pair separated 
by N transverse links will be connected by exactly N glu- 
ons — one on each links in between 0. The reason for 
doing this approximation is that it allows to illustrate the 
confinement mechanism on the LF more clearly. Further 
below wc will show results that do not make use of this 
approximation. 

The physical meaning of the one gluon per link approx- 
imation is as follows: In the limit of large gluon masses 
(=large lattice spacings) a QQ pair is connected by a 
static (no fluctuations) chain of gluons. As the gluon 
mass gets smaller, the string starts fluctuating. Fluctua- 
tions where the string moves longitudinally with respect 



^Note that even though this has been shown only for 
dimensional gauge theories in Ref. , the result can be used 
here since QCD on a transverse lattice is formally equivalent 
to a large number of coupled 14-1 dimensional field theories! 

^Without this approximation, this would be the minimal 
configuration, i.e. any physical state subject to the color sin- 
glet constraint at each site will have at least this number of 
gluons. 
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to the straight line connecting the QQ pair correspond to 
excitations which do not change the number of hnk fields. 
Those are included in the above approximation. Fluctu- 
ations which are excluded in the above approximation 
are the ones where the string gets deformed transversely 
so strongly that it winds forth and back in the trans- 
verse direction. If one wants an exact solution to QCD, 
all fluctuations must be included. However, as a first 
step beyond the (static !) heavy gluon limit, it makes 
sense to include only the "motion" of the gluons first 
but not creation and annihilation of additional gluons. 
Note that if one were to work close to the continuum 
limit then such an approximation would not make sense. 
However, since we consider the Vs as blocked variables, 
which correspond to rather rigid degrees of freedom, it is 
not completely unreasonable to assume that pair creation 
is suppressed in ground state configurations. 

The general formalism for calculating the (rest-frame) 
QQ potential within the light-front framework has been 
discussed in detail in Ref. and the reader in strongly 
encouraged to consult these works for details. In the fol- 
lowing we will only present the explicit integral equations 
that one obtains for the coupled QQ — glue system. If 
the quark and the anti-quark are on the same site, with 
longitudinal separation x~ ,the calculation is trivial since 
the " integral equation" collapses into one single equation. 



P- 



(3.2) 



where = g^Cpli^ and Cp = {Nl - l)/2Nc. g car- 
ries dimensions of mass. In the naive continuum limit, is 
related to the coupling constant go in the Lagrangian of 
QCD2+1 via g^ (X go/a, where a is the transverse lattice 
spacing, which is why the coupling has this unusual di- 
mension. As has been shown in Rcf. ||l3|, x~ is related 
to the longitudinal separation of the QQ pair in its rest 
frame xl via 



XL 



X V 



(3.3) 



where v'^ is the velocity vector of the QQ pair. Similarly, 
is related to the potential interaction energy between 
the quark and the anti-quark (again in their rest frame) 
through the relation 



V{xl,xa 



P- 



i.e. for zero transverse separation one obtains 
V{xl,0)^G'^\xl\. 



(3.4) 



(3.5) 



As soon as the QQ pair is separated by at least one trans- 
verse site, the result is less trivial. For example, when the 
transverse separation is one lattice spacing, one obtains 
an integral equation for the wave function of the (one !) 
gluon connecting the QQ pair 



P~xP{k^ 



k+ 
2^ 



TO 

2k+ 



V'(fc+) 



(3.6) 



G2 / dq^ 



{q+ + k+) U{k+) - i/;(g+)e»(''^+-9+)^"/2 



G2 



(9^ 



2y/k+^ {q+ ^ k+y 
) [V(fc+) - V(9+)e'('^"''^^"="/^ 



2y/k+^ {q+ - k+Y 



where the two interaction terms arise from the Coulomb 
coupling of the gluon t o th e quark and anti-quark respec- 
tively. Eqs.(p^) and (3.4) hold also here. The exponen- 
tial "form -factors" in the interaction terms arise since 
quark and anti-quark are displaced in the longitudinal 
direction. 

In the continuum (i.e. when one solves this integral 
equation exactly) the resulting V(xltX±) is independent 
of the velocity w"*". Note that this is not the case when 
one uses DLCQ to solve the integral equation. This point 
will be discussed below. 

For a separation of two sites the integral equation for 
the two gluons between the charges reads 



p-^{kt,k+) = 



G' dq^ 
Jo 



G'l dq^ 
la 



2k+ 2k+ J ^ ^ 2 
{q+ + k+)\tlj{k+,k+) - V'(q+, fc+)e'('=i'-'?^ 



2v+^ 



/cV q 



+ 



1 y \f 



{q+ + k+) ^{k+,k+)-^{k+,q+)e 



-kty 



G- 



dq+- 



2^fkU- 
kt) {k+ + 



-kt) 



i^k+k+q+{kt + k+-q+) 

[^{k+,k+)-^iq+,k+)] 



{q+ - k+ 



G^ 



Ay/kJ^; 



-.^{k+,k+) 



(3.7) 



The three interaction terms in Eq.(3.7) arise from the 
interaction of the first gluon with the quark, the second 
gluon with the anti-quark and the interaction between 
the two gluons, respectively. 

The generalization of these expressions to more than 2 
links is straightforward but the resulting expressions are 
very lengthy and will be omitted here. 

Note that we used, the "Coulomb trick" ||25|l in 



Eqs.(3^) and (3/7) by adding and subtracting analyti- 
cally a term in the interaction. This results in an in- 
teraction that vanishes for constant wave functions — 
which is close to the actual shape — and thus numerical 
convergence is improved considerably: 



dy ] {x + y){2-x-y) 
{x- yf \ ^/xy{l - x){l - y) 



V'(y) - 4-0(2;) 
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dy ] {x + y){2-x- y) 



Jo (x - y f \ ^/xy{l - x){l - y) 
mplx) Atplx) 

as well as 



[V'(y) - ■^[x]] 



(3.8) 



dy 



- y) 

- dy 



{x + y) 



{x-y) 



xy 

{x + y) 



/xy 



(3.9) 



Even though these integral equations will be solved be- 
low in their fully relativistic form, it is very instructive to 
consider various approximations thereof — particularly 
the limit iv? — > oo. Making a non-relativistic expansion 
around the minimum of the kinetic terms one finds for 
the (rest frame) potential energy for an arbitrary config- 
uration with gluon (rest frame) positions at xi, ...,a;„j^ 



V 



XL 



-xi\ 



XL I 

2 ' 



+ J2 



Xi+l\ 



(3.10) 



Clearly, this expression is minimized for 

<a;i<...<a::„i < ^ with minimum value 



XL 

2 



V{xl,x±) = n±m + G'^ — \xl\- 



(3.11) 



For large m this is the QQ potential in the rest frame 
of the pair, which thus exhibits linear confinement 
Since the lattice spacing has not yet been fixed so far, we 
are free to choose 



2m 



(3.12) 



which renders the longitudinal and the transverse string 
tension equal to each other. Nevertheless, V in Eq.( 3.11 ) 
is still not rotationally invariant: for example, for xl = 
x± = X one finds V{x,x) = 2V{0,x). For a rotation- 
ally invariant linear potential the result would have been 
V{x,x) = V2V{0,x). 

Such a result is familiar from the strong coupling limit 
in Euclidean or Hamiltonian lattice QCD, where one ob- 
tains exactly the same square shaped equipotential lines. 
Of course, once one no longer restricts oneself to the 
strong coupling limit, the quantum mechanical fluctu- 
ations tend to restore rotationally invariance. 

This is also what happens here. To see this, let us 
consider a QQ pair separated by one transverse lattice 
unit. When the (rest frame) positions of the quark and 
the anti-quark are both at xl = then the gluon in 
between them experiences a potential energy equal to 



Valignedixi) = G^TT 



(3.13) 



In contrast, when the positions of the quark and the anti- 
quark are at ±^ ^ the potential energy which the 
gluon in between them experiences is 



^displaced {xi ) 



G^TT 



Xl_ 

' 2 



Xl 



xl_ 
2 



Xl 




(3.14) 



The violation of rotationally invariance in the large m 
limit manifests itself through the fact that the minimum 
of the potential (w.r.t. Xi) for the displaced case is too 
high compared to the aligned case. 
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FIG. 2. Potential energy of a gluon between an ex- 
ternal QQ pair separated by one link in the (fictitious) 
non-relativistic limit, a.) for Q and Q at the same longitu- 
dinal position, b.) for Q and Q displaced in the longitudinal 
direction hy xl = 2/G. 

The crucial point is ( see Fig.^) t hat t he valley of the 
potential in the displaced case [Eq.(3.14)] is wide r tha n 
the valley of the potential in the aligned case [Eq.( 3.13 )]. 
Therefore, corrections due to the quantum mechani- 
cal zero point energy tend to increase the energy more 
strongly in the aligned case than in the displaced case. 
In a sense, quantum mechanical corrections work in the 
right direction to help restore rotationally invariance. 

Even though this simple quantum mechanical argu- 
ment is very useful in order to understand the physics 
of the restoration of rotationally invariance, we will not 
attempt to make it quantitative since we will now make 
a fully relativistic calculation, where we actually solve 
light-front integral equations for these systems, Eqs. 
( |3.6D , (3.7) and the generalization to more sites. 

As we have seen above, the mass term plays an im- 
portant role in setting the scale for the transverse string 
tension. In fact, given the string tension a, one obtains 
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a = <j/m for large m? . Thus, if one wants to be close to 
the continuum limit, one should try to make rn? as small 
as possible. It turns out that below = 0, the spec- 
trum becomes tachyonic i.e. m? = is the smallest 
meaningful value. As we will see below, the string tension 
in lattice units (i.e. also the lattice spacing in physical 
units) remains finite at m? — 0. Since m? = is thus 
the closest we can get to the continuum limit within our 
approximation, we will focus only on this value. 

In order to fix the scales, we first determine the trans- 
verse string tension. There are two possibilities to de- 
termine this observable. One is to consider the QQ po- 
tential for large separations of the pair. The other is 
to consider a torus geometry, where the energy of glue- 
ball states which "wrap around" the torus is equal to the 
string tension times the circumference of the circle [ p4[ . 
Even though the latter method cannot be used to calcu- 
late the QQ potential or to calculate the string tension 
in any direction other than the transverse direction, it 
turns out to converge faster than calculations of the QQ 
potential. Therefore, " wrap-around glueballs" were used 
to fix the string tension in this work. 

In the calculation of the glueball masses, we used 
DLCQ with anti-periodic boundary conditions, be- 
cause the convergence in the DLCQ parameter K is faster 
than with periodic boundary conditions. Finally, the 
masses of "wrap-around glueballs" were calculated for 
a fixed size of the periodic lattice as a function of the 
DLCQ parameter K and the results were extrapolated 
to if ^ oo (Fig.|). 
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FIG. 3. Ground state energy of winding modes for several 
transverse sizes n± of the lattice as a function of the inverse 
DLCQ parameter K. The shaded bands correspond to the 
isT — > cxD extrapolated results for the ground state energies. 
The width of the bands reflects systematic uncertainties in 
the extrapolation. 

The transverse string tension (in lattice) was then ex- 
tracted by considering the resulting extrapolated masses 
of these states as a function of the number of lattice 
spacings. Note that, as the nearly equidistant spacing 
in FigJ^ indicates, already for only a few lattice spacings, 
the mass of these states depends almost exactly linearly 
on the number of transverse lattice spacings, which al- 
lows one the easily extract the transverse string tension 
from this data. The fitted result is 



M{n) 
G 



ncr_La« n (1.30 ±0.03) , 



(3.15) 



where we introduced the transverse string tension and the 
transverse lattice spacing. If one denotes the longitudinal 
string tension by ctl, one finds (|3.5|) 



ctl = G — 



(3.16) 



and if one furthermore imposes ct^ = ctj^ = cr one finds 
for the transverse lattice spacing in physical units 



def lim„- 



M{n) 



^(0.83 ±0.02) (3.17) 



Having related the longitudinal and transverse scales, we 
are now in a position to explore the QQ potential. Note 
that even though the "wrap-around glueballs" are quite 
useful in the determination of the purely transverse string 
tension, they cannot be used to determine the string ten- 
sion for any direction other than the transverse direction 
and they also cannot be used to determine the QQ poten- 
tial for finite separations. The QQ potential will there- 
fore be determined by using the fixed charges formalism 
outlined above. Nevertheless, since — as we will see be- 
low — the determination of the string tension from the 
QQ potential through the fixed charges formalism is less 
accurate than the method using wrap-around glueballs, 
we will still use the latter method to fix the scales in the 
following calculation. 

In the calculation of the QQ potential we proceeded 
as follows: For a given longitudinal separation x'^ and a 
given velocity u+, the DLCQ Hamiltonian (for a given 
number of transverse sites) was constructed with a cutoff 
on the total longitudinal momentum of the gluon string^ 



^In contrast to DLCQ calculations of glueball masses, mo- 
mentum is not conserved here since arbitrary momenta can 
be transferred to and from the external charges. 



5 



Since the momenta are discrete, provides an infrared 
cutoff in the following sense: the peak of the wave func- 
tion, will be at momenta of order fc+ w u+G. I.e. only 
if v^G is large in integer units, one will not be affected 
by the cutoff. Thus one must perform a careful extrapo- 
lation^, were one sends both the UV cutoff, but also the 
velocity to infinity. The actual numerical work was 
based on DLCQ with anti-periodic boundary conditions 
on the link fields. The ground state energies and wave 
functions for the resulting DLCQ Hamiltonians were de- 
termined using a Lanzcos algorithm p6[ |. The resulting 
QQ potential is shown if Fig.|[ 
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lines in Figjj. The residual anisotropy is mostly d ue t o 
the use of an oversimplified effective potential [Eq.( ^.4[ )] 
but also due to the (ad hoc) suppression of higher Fock 
components (see also the results below obtained without 
Fock space truncation). Nevertheless, the restoration of 
rotational invariance is quite impressive as a contour plot 
for the same data as in Fig.^ shows (Fig.||). 




-3-2-10 1 2 3 




FIG. 4. QQ-potential versus r, where = x\+ x\, as ex- 
tracted from the light front Hamiltonian in the approximation 
where only one gluon per link is allowed. The error bars reflect 
the uncertainties arising from the extrapolation to the contin- 
uum limit in the DLCQ calculation. The transverse scale was 
fixed from the masses of "wrap-around glueballs". The fact 
that not all points lie on the same smooth curve reflects the 
residual anisotropy due to the Fock space truncation and the 
use of an oversimplified efi'ective potential. For comparison, 
the two dashed lines give the range of values one would ob- 
tain in the large m limit, 1/ oc |a;L| -I- \xi_\. Dots, crosses, 
diamonds, squares and circles correspond to n±_ = 0, 1, 2, 3, 4 
respectively. 

Even though not all points in Fig.^ lie on the same 
smooth curve, there is still a significant improvement 
compared to the large m limit (3.11), where values for 
V{xl, a^j.) would fill the whole area within the two dashed 



*This must be done very carefully because the peak of the 
wave function moves as one changes v'^ . 
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FIG. 5. a.) Contour plot for the QQ-potential ver- 
sus x± and XL in the approximation where we keep only 
one gluon per link. For comparison, the contour plot for 
the QQ-potential in the large m (strong coupling) limit 
{V oc \xl\ + \x±\) is shown in b.). 

The lines of constant potential turn out to be almost 
circles. The slight anisotropy of the data in Fig.|| results 
in a slight ellipsoidal distortion in Fig.^]. Note that the 
ellipsoidal shape in Fig|| does not mean that the string 
tension in the longitudinal and transverse directions are 
different! In fact they are the same by construction, i.e. 
by the choice for the lattice unit. What is going on here 
is that, while V{xl,0) — G'^^\xl\ in the approximation 
chosen in this work, the potential in the transverse direc- 
tion satisfies V{0,xj_) ^±zl^ G^^\x±\ +c, which has the 
same string tension in longitudinal and transverse direc- 
tion but which results in a equipotential lines that look 
like ellipses. 



IV. THE MOMENTUM CARRIED BY THE GLUE 

Even though the calculations presented in this paper 
are still very crude, the approximate rotational invari- 
ance of the QQ potential is very encouraging so that we 
proceed to investigate other physical observables. We 
picked the gluon distribution in these QQ systems since 
one of the main motivations to study light-front QCD in 
the first place is the direct access to parton distributions 
measured in deep inelastic scattering. In particular, we 
will focus on the (light-front) momentum carried by the 
gluon component. 

Before we proceeds, we should caution the reader that 
the results must be interpreted with care: Since we work 
with finite and large lattice spacing, the "gluons" are 
not point-like "current gluons" and thus the distribution 
functions should be considered as distribution functions 
of some kind of "constituent gluons". Nevertheless, the 
results might be helpful in obtaining some intuitive in- 
sight about the parton structure at low Q^. 

Obviously, when the external charges are sitting on the 
same site, there are no gluons within the approximation 
considered in this paper and thus the gluons carry zero 
momentum for such a configuration. The other extreme 
is when the charges are separated only transversely, in 
which case one can derive an exact result: 





Consider Eq.(3.6) for x 



^/;(fc^ 



(4.1) 



_^g2 n^+ (g++fc+)[^(fc+)-vxg+)] 



2^k+q+ {q+ ^ k+f 



n + (g+ + fc+)[v^(fc+)-V'(g+)] 

Jo 2y/k+^ {q+ - k+f 

On the one hand (see Appendix A) one knows that 



P- 



(4.2) 



where {V) is the expectation value for potential (rest- 
frame) energy in this configuration. On the other hand, 
from the Feynman-Hellma n th eorem for variations of 
with respect to in Eq.(4.1) one finds 



-r^p- = ^ / dk+ \Mk^ 

dv+ Jo ' 



■k+. 



(4.3) 



Combining these two results, one thus finds the remark- 
ably simple result (valid for QQ pairs that are separated 
only transversely) 



dk+ \^{k+)\^ k+ = {V)v+ (4.4) 



which states that the momentum carried by the gluons is 
proportional to the energy in the gluon field. Applying 
the same reasoning to Eq.(3.7) (again for = 0), one 
obtains similarly 



P. 



glue 



dk+dk+\i,{k+,k+f {kl 



{V)v+, 
(4.5) 



which has the same simple interpretation as Eq.( [4.4[ ). 
The generalization to more than 2 links is straightfor- 
ward and obvious: 



3+ — 



dk+...dk+ \^P{k+, k+f {k+ + ... + k + ) 



= {V)v^ 



(4.6) 



Note that the similar sum rules have been derived else- 
where and they are not limited to the Fock space 
truncated version. 

In the general case, where both xj^ and x± are nonzero, 
we were not able to find a generalization of Eq.(4.6). 
However, we can easily use the numerically obtained 
states from the calculation of the ground state energy 
and just "measure" the momentum carried by the glue. 
The result is shown in Fig.H. 



7 



"3) 



" 2 

o 1 



0.0 



0.5 



1.0 



1.5 



2.0 



2.5 



3.0 



3.5 



XlG 



FIG. 6. Momentum carried by the gluons for external 
charges that are separated by n± = 1,...,4 lattice spac- 
ings as a function of the longitudinal (rest-frame) separa- 
tion of the external charges. The results are obtained in the 
one-gluon-per-link approximation. For n± = the gluons 
carry no momentum in this approximation. 

First of all, one notices that the momentum carried 
by the gluons increases as the transverse separation in- 
creases. This result is intuitively obvious and happens 
for XL = in a way that is consistent with Eq.(|4.(j). A 
much less obvious result is that, for fixed nj_, the mo- 
mentum carried by the gluons decreases for increasing 
longitudinal separation. This result is at first counter- 
intuitive because the energy in the gluon field increases 
when one increases xl for fixed n±. However, this ap- 
parent paradox gets resolved when one recalls that only 
transverse components of the electric field contribute to 
the momentum of the state in the infinite momentum 
frame |2^. Assuming (for simplicity) that the QQ pair 
is connected by a straight gluon string, it is obvious that 
the transverse component decreases as the string gets ro- 
tated into the longitudinal direction. This also explains 
why the momentum carried by the gluons depends not 
only on the separation of the QQ pair but also on its ori- 
entation. As an idealized case, it is instructive to assume 
that the QQ pair is connected by a straight string, i.e. 
the electric field lines are parallel to the QQ-axis (which 
has an angle with the longitudinal direction that we will 
denote by 0. If is the field strength in the string then 
\E±\ = \E\ sin 6. In Ref. it has been shown that the 
momentum carried by the gauge bosons is proportional 
to the volume integral of i?^, i.e. one would expect 



^aiuel^^ oc r sir? Q = x±^ sin Q 



(4.7) 



which qualitatively explains the decrease of P^i^^^ for fixed 
x±_ with increasing xl- 

The qualitative conclusion that one can take from these 
results is that for heavy QQ systems, such as J /'^ or 
Upsilon mesons, it is quite natural that the gluons carry 
a fraction 



^ _ {V) 



total 



M, 



QQ 



(4.8) 



of a hadron's momentum, where c is a " geometry factor" , 
which is of the order 1. 



V. INCLUDING HIGHER FOCK COMPONENTS 



Even though the results in Section |[I| are already 
amazingly close to being rotationally invariant, they have 
a serious flaw: they were obtained after an ad hoc trun- 
cation of the Fock space to the minimal component (< 1 
quanta per link). While this provided us with a very 
intuitive and numerically acceptable "quantum mechan- 
ical" approach to the gluon distribution, it gives rise to 
results that are still very close to the strong coupling 
limit: in fact, even though the fields may fluctuate longi- 
tudinally (which distinguishes the above valence approx- 
imation from the strict strong coupling limit) QED and 
QCD still give identical static potentials at this level. 
One might thus be sceptical about the results obtained 
— especially about the linearly rising QQ potential. In 
this section, we will provide evidence that the results ob- 
tained for the QQ potential do not depend very much on 
the valence approximation. 

The main difference between the results in this section 
and the results in the previous section is that we allowed 
for one additional link-antilink pair in the string state. 
Computer storage and the algorithm that we were using 
(Lanzcos, where we stored all nonzero matrix elements 
and their addresses in order to keep the code fast) did 
not allow us to go beyond one additional pair while at 
the same time being able to extrapolate to the continuum 
limit for two or more links. However, our calculations 
with one additional pair showed that this component of 
the wave function is already very strongly suppressed 
compared to the valence component so that we do not 
expect significant changes of the result by going to even 
higher Fock components. Also, we should point out that 
with just one additional pair there is a difference between 
QED and QCD. 

The effective potential was taken to be the same as 
the one we used for the valence calculation, namely just 
a quadratic term and we did the calculations with the 
"mass" at the 1** order critical point. We should empha- 
size that adding higher order terms to the effective poten- 
tial only affects the calculation through higher Fock com- 
ponents. Therefore (unless the coefficients of the higher 
order terms are taken to be very large to invalidate the 
perturbative argument) we expect the QQ potential to 
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be rather insensitive to such higher order terms. This is 
in sharp contrast to gluebaU calculations where even the 
ground state spectrum depends strongly on the higher 
order terms j24]. The rest of the calculational procedure 
strongly resembled the calculation in the previous sec- 
tion. Because of limitations of our algorithm, we had to 
restrict ourselves to 3 or less lattice sites. 

The most important result was that higher Fock com- 
ponents in the ground state of the string had a norm of 
typically only a few per cent, with the resulting energy 
shifts slightly larger but still of the same order of mag- 
nitude. In fact, this justifies the truncation to only one 
pair. The results for the QQ potential and the momen- 
tum carried by the glue are shown in Figs. (0) and (||) 
respectively. 
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FIG. 7. Contour plot for the QQ-potential versus x± and 
XL the approximation where we allow for up to one addi- 
tional link-antilink pair. 

The first thing that one notes when one compares Fig. 
^ and Fig. ^ is an increase in the size of the error bars. 
This is because we had to limit ourselves to smaller values 
of the longitudinal momentum and the systematic uncer- 
tainties from the extrapolation increased. Another obvi- 
ous observation is that by including an additional pair the 
result gets even closer to being rotationally invariant. We 
find that very encouraging for future calculations. Apart 
from that, the only difference between the calculations 
with and without the extra pair is a renormalization of 
the string tension in lattice units (aj^ = 1.14(t~^/^ com- 
pared to a_L{valence) — 1.04(t~^/^). 

The momentum carried by the glue no longer has to 



vanish for a QQ pair that is oriented longitudinally (Fig. 
^). Nevertheless, the numerical result is very small, 
which reflects the small admixture from higher Fock com- 
ponents. In fact, similar to the potential, there is only 
very little change for the momentum carried by the glue 
when one compares calculations with (Fig. ^ and with- 
out (Fig. ^) link-antilink pairs. 
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FIG. 8. Momentum carried by the gluons for external 
charges that are separated by n± = 0, ...,3 lattice spac- 
ings as a function of the longitudinal (rest-frame) separa- 
tion of the external charges. The results are obtained in the 
one-gluon-per-Unk approximation plus one additional pair. 

Overall, we found only a small change in the results for 
the ground state of the string as we allowed for higher 
Fock components. We find this a very encouraging result 
since it may allow one to construct models for quarkonia 
that have a similar truncation in the Fock space. 



VI. SUMMARY AND OUTLOOK 

We have presented analytical and numerical calcula- 
tions for the rest-frame QQ potential in QCD on a trans- 
verse lattice using light-front quantization. The rest- 
frame QQ potential was obtained by considering the in- 
variant mass of a heavy QQ pair that moves with con- 
stant velocity and fixed separation. 

First we showed analytically, that one obtains a lin- 
ear potential (both in the longitudinal as well as in the 
transverse direction) in the limit of large transverse lat- 
tice spacings. This limit is similar to the strong coupling 
limit in Hamiltonian or Euclidean lattice QCD. The con- 
finement mechanism for QCD in LF gauge in this limit 
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depends on the orientation of the external charges: In 
the longitudinal direction, confinement arises from the 
instantaneous interaction which appears when the un- 
physical component of the gluon field is eliminated by 
solving the LF gauge analog of the Poisson equation. 
This component of the gluon field is also unphysical in 
the sense that it carries no momentum in the infinite 
momentum frame. Thus, for a purely longitudinally sep- 
arated QQ pair the gluon field carries no momentum at 
all in the limit of large transverse lattice spacings — 
even though there is potential energy in the gluon field. 
For transversely separated charges, a completely differ- 
ent confinement mechanism is at work for large trans- 
verse lattice spacings: gauge invariance demands that 
the QQ pair is connected by a string of gauge-link field.^ 
The "strong coupling limit" in Hamiltonian QCD cor- 
responds on a transverse lattice to the limit, where the 
Lagrangian for the gluon link contains a large mass and 
thus the ground state energy of a gluon configuration is 
obtained by counting links. Together with the instanta- 
neous interaction between quarks and adjacent link fields 
as well as among adjacent link fields, one thus obtains a 
square shaped linear potential — a result that is familiar 
from the strong coupling limit. Note that the link fields 
carry momentum in the infinite momentum frame, i.e. 
for transversely separated QQ pairs, the gluons do carry 
a sizeable fraction of the total momentum. 

Even though such an extreme asymmetry — no mo- 
mentum carried by the gluons for purely longitudinal sep- 
arations — is the result of the strong coupling limit, it is 
actually quite natural and physical that there is such an 
asymmetry depending on the orientation. This effect oc- 
curs already in QED and arises from the fact that the 
component of the (color-) electric which is transverse to 
the boost direction transforms differently from the com- 
ponent parallel to it. In particular, the parallel compo- 
nent does not contribute to the Pointing vector in the 
infinite momentum frame. For this reason, and because 
of the difference in the field distribution between QED 
and QCD (dipole versus string) we expect a stronger ori- 
entational dependence of the gluon momentum in QCD 
than in QED. 

In the next step we performed numerical calculations 
of the QQ potential. First we showed results where we 
truncated the Fock space to not more than one quantum 
of the link field for each transverse link. Otherwise, the 
link fields were treated fully dynamical and were allowed 
to move freely in the longitudinal direction. There was 
only a first order transition as a function of the link field 
mass, i.e. while the transverse lattice spacing in physical 
units decreased with the mass term, it did not approach 



zero at the critical point. Nevertheless, the numerically 
obtained QQ potential at the critical point was almost 
rotationally invariant. The reason for the first order crit- 
ical point was most likely our oversimplified choice of the 
effective link field potential Veff{U). 

Even though the numerical precision was only limited, 
we then showed that the qualitative results change only 
little when we allowed for excited Fock space compo- 
nents. This is very important, since only upon inclusion 
of higher Fock components is there a difference between 
QED and QCD. The unimportance of higher Fock com- 
ponents in the ground state of the string is presumably 
partly responsible for the approximate rotational invari- 
ance of our results. The reason is that in general one 
would not expect rotational invariance on a coarse lat- 
tice unless one includes the full effective potential for the 
link fields — from which we kept only the quadratic term. 
Those terms that we omitted, contribute only to matrix 
elements that include higher Fock components. There- 
fore, if higher Fock components are not important for 
the QQ potential then the higher order terms in 14// (C/) 
cannot be important either.^ For tests of rotational in- 
variance that are sensitive to the higher order terms in 
Veff{U) one must therefore study other observables, such 
as glueball spectra p^ . 

The numerical results, both with and without trunca- 
tion of the Fock space, confirmed the strong dependence 
of the momentum carried by the glue on the orientation 
of the external pair. This might have observable conse- 
quences for heavy quarkonia with nonzero orbital angular 
momentum. 

There are many extensions of this work that one could 
think of, such as studying excited states of the string or 
using Monte Carlo algorithms for obtaining the eigen- 
states and eigenvalues of the LF Hamiltonian. Excited 
states would be interesting both for theoretical but also 
for phenomenological reasons. On the theoretical side, 
since higher Fock components play a stronger role in the 
excited states of the string, they would allow probing the 
relevance of the higher order terms in 14// (t/) for rota- 
tional invariance and could thus be used to help fix those 
constants. From the phenomenological point of view, the 
excited states of the string are interesting because of their 
connection to hybrid states. Employing Monte Carlo al- 
gorithms jsot when studying the transverse lattice might 
be a useful option since that would allow one to include 
many more Fock components than otherwise. Repeat- 
ing the calculations in 3 -I- 1 dimensions is an obvious 
extension of this work . One possibility in this direc- 
tion, which we are currently exploring, is to repeat the 



^Even if one is not dogmatic about gauge invariance, in- 
frared divergences, that occur otherwise, practically enforce 
this condition. 



''Note that this argument is not fully complete since we have 
not verified that the higher Fock components are unimportant 
for the ground state of the string after the higher order terms 
are included in Ve//(I7). 
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minimal Fock space truncation and to fit the plaquette 
interaction term such that approximate rotational invari- 
ance in all 3 spatial directions is achieved. Then one can 
use this gluonic interaction to investigate the physics of 
the ground and excited states of heavy quarkonia. 
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VII. APPENDIX: HOW TO MEASURE THE 
REST-FRAME QQ POTENTIAL IN A 
LIGHT-FRONT CALCULATION 

In this Appendix^ a scheme is developed which allows 
one to extract the QQ potential, i.e. the quantity which 
corresponds to the potential between two infinitely heavy 
quarks in a rest frame, from a light-front calculation.^ 
There are several reasons to study this observable in the 
light-front (LF) framework 

• The LF formalism lacks manifest rotationally in- 
variance. Therefore, if one starts with a wrong LF 
Hamiltonian for QCD, the result V{R) depends on 
the orientation of R with respect to the 3-axis. ^ A 
measurement of V{R) thus provides a direct probe 
of rotationally invariance in a physical observable. 

• This sensitivity of V{R) to rotationally symme- 
try can then be exploited in the renormalization 
procedure to help determine non-covariant counter 
terms. 

• And, most importantly, V{R) in QCD is very well 
known over a large range of distances from the spec- 
troscopy of heavy QQ mesons as well as from non- 
perturbative Euclidean lattice calculations (at least 
in the absence of dynamical quarks — but it is easy 
to make the same approximation in a LF frame- 
work). 



''See also Ref. ||T|l. 

*Even though the body of this paper in on 2 + 1 dimensional 
QCD, we keep the discussion in this section independent of 
the number of transverse space-time dimensions. 



We use the notation A 
{A\A^). 



A^ = (A« ± ^^)/^/2, A^ 



FIG. 9. World lines for two charges with longitudinal sep- 
aration Aa;^ in the rest frame. 

Before we embark on deriving the effective LF Hamilto- 
nian for two infinitely heavy sources, it is instructive to 
understand physically what it means to have two fixed 
sources "at rest "p^ from the LF point-of-view. As should 
be clear from Fig^, fixed charges in a conventional frame 
correspond to charges that move with constant velocity 

on the LF ~ ~ 1/a/2 in the example in Fig.^. 
Furthermore, if the longitudinal separation is Ax^ in the 
rest frame, the charges have fixed separation Ax^ = 
^/2Ax^ in the longitudinal LF direction. In the more 
general case, where the charges are moving with con- 
stant four velocity w'', where v± = in the rest frame, 
one obtains Ax" = Ax^ /v~^. 

The transverse separation is the same on the LF as 
it is in a rest frame description. Therefore, in order to 
understand the LF-physics of charges that are fixed in a 
conventional frame with separation R = {R-^, R^,R^), we 
must first understand how to describe a "dumbbell," with 
ends separated by (Aa;\ Ax^, Ax") = i^^ 
that moves with constant velocity w^. 



A. One Heavy Quark on the LF 

A pair of sources, moving both with the same constant 
velocity can also be interpreted (and treated) as one 
extended source moving with constant velocity v'^ (for 
simplicity, we will keep v± — 0). This is reminiscent of 
heavy quarks and thus, as a warmup exercise, it is very 
instructive to consider one (point-like) heavy quark on 
the LF first (see also Refs. |l|j2§). 

For simplicity, we will first take the heavy quark limit 
for the canonical Hamiltonian, which can be written in 
the form 



^"Here we mean "at rest" for a conventional observer, i.e. in 
an equal time frame. 
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/l/f2 + P 



(7.1) 



where _B represents the hadron, b is the heavy quark, P^j^ 
contains the interactions between heavy (6) and hght de- 
grees of freedom and contains all terms involving 
light degrees of freedom only. The heavy quark limit is 
obtained by making an expansion in inverse powers of 
the 6-quark mass. For this purpose we write 



Pt = Pb - Pt = Mbv^ 



Pt 



(7.2) 



where is defined to be the sum of the longitudinal 
(LF-) momenta of all light degrees of freedom. For the 



(total) LF-energy we write on the l.h.s. of Eq.(7.1) 



P„ = Mbv- = 



Mb 
2v+ 



Mb + SE 
2v+ ' 



(7.3) 



where SE = Mb — Mb is t he " bindi ng e nergy" of the 
hadron. After inserting Eqs.(7.2) and (7.3) into Eq.(7.1) 
and expanding one obtains 



Mb + 6E 
2v+ 



2{Mbv+-p+) 



P, 



HL 



P, 



LL 



Ml 



2 {MbV+ + SEv+ - p+ 

Mb 5E 
2^ ^ 2^ 



P, 



HL 



P, 



LL 



2w+2 



O 



1 



P, 



HL 



P; 



LL- 



gluon exchange interactions with light quarks or glu- 
ons, which are, respectively, Vgq oc (p^ — Pq^)~^ and 

Vqg cx (pg + p'g^ ) {p'^ — p'g^ ) ^ and remain unchanged in 
the limit Mb — > oo. Terms which involve instantaneous 
gluon exchange and are off-diagonal in the brown muck 
Fock space behave in the same way. 

The quark gluon vertex simplifies considerably. For 
finite quark mass one has for the matrix element for the 
emission of a gluon with momentum k, polarization i and 
color a between quarks of momentum pi and p2 



P 



QQg 



, 2 A;! _ 



,- g_LPi_L + iMb 
Pt 
(7.6) 



where spinors as well as creation/destruction operators 
have been omitted for simplicity. In the heavy quark 
limit [note l/p+ - l/p+ = 0{M-^)] the spin dependent 
terms drop out and one finds in the heavy quark limit 



k+' 



(7.7) 



(7.4) 



The spin of the heavy quark thus decouples completely, 
giving rise to the well known SU{2Nf) symmetry in 
heavy quark systems. 



B. Two Heavy Sources 



Note that we have assumed that the transverse momen- 
tum of the heavy quark is small compared to its mass, 
which is justified in a frame where the transverse velocity 
of the heavy hadron vanishes. The term proportion al to 
Mb cancels between the l.h.s. and the r.h.s. of Eq.(7.4) 
and we are left with 



5E 



Pl 
2w+2 



- Oil/Mb) + Pj 



HL 



P, 



LL- 



(7.5) 



The brown muck Hamiltonian P^^ is the same as for 
light-light systems and will not be discussed here. The 
interaction term between the heavy quark and the brown 
muck (P^l) more tedious but straightforward. For 
example, heavy quark pair creation terms (via instanta- 
neous gluons) are proportional to l/(p^^ +P;^)^ oc l/M^ 
and can thus be neglected. Similarly, pair creation of 
heavy quarks from virtual gluons is also suppressed by 
at least one power of Mb- This also justifies our omis- 
sion of states containing more than one heavy quark 
from the start. Other terms that vanish in Pff^ in- 
clude interactions that involve instantaneous exchanges 
of heavy quarks, which are typically proportional to the 
inverse p"*" of the exchanged quark and thus of the or- 
der 0{Mf^^). Up to this point, all interaction terms 
that we have considered vanish in the heavy quark limit. 
The more interesting ones are of course those terms 
which survive. The simplest ones are the instantaneous 



As we discussed above, two heavy sources at fixed sep- 
aration can formally be treated as one extended heavy 
source |^. Therefore, the LF Hamiltonian for two sources 
is the same as for one source with two minor modifica- 
tions: 

• All vertices involving a heavy source get modified 
according to the rule ( O stands for any operator 
acting on the brown muck) 

{x^'t"x X O} + /i.e. ^ (7.8) 

{ [xh'T'XQFRiq) - X^QT'^XQ'FB.i-q)] x o} + h.c, 



where x stands for the operator acting on the color 
degrees of freedom of the heavy quark/ antiquark. 
The "form factor" 



FR{q) = exp 



1 fq+R^ 

2 [ v+ 



R±qA 



(7.9) 



arises from acting with the (kinematic !) dis- 
placement operator on the position of the heavy 



Just think of a dumbbell. 
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quark/antiquark (shifting it from x =0, x±^ = Oj^ 
to x~ = ±i?^/2w+, x_L = ±R_l/2) and q is the net 
momentum transferred to the brown muck. This 
rule holds irrespective of the number of gluons in- 
volved in this process. Note that 0{brownmuck) 
is the same for one or two heavy sources. 

• There is a static potential between the two heavy 
quarks. In the continuum, the canonical Hamil- 
tonian yields 

Phh = g^X^qT-XQX^T-XQ' ^^^KR^mv^. In 
general, there will be a more complicated depen- 
dence on R which has to be determined by demand- 
ing self-consistency. For example, singularities aris- 
ing from exchange of gluons with low g+ between 
the two sources should cancel (nonperturbatively) 
with the IR behavior of the instantaneous potential 



HH 



Using Eq.(7.5), V{R) can thus be extracted as follows: 



1. For a given R and write down the effec- 
tive LF-Hamiltonian for the heavy pair interacting 
with the brown muck (including the form factors 
Eq.( [7.9| ) and including all the counter terms and 
counter term functions which would also appear in 
a "heavy-light" system). 



2. The lowest eigenvalue 5E^^^ from Eq.(7.5), i.e. the 
QCD-ground state in the presence of the two heavy 
sources, is then equal to V{R). 

The heavy quark potential thus calculated is equivalent 
to the potential which a lattice theorist would extract 
from an asymmetric rectangular Wilson loop. Since there 
is plenty of "quenched" lattice data around, it would 
make sense to omit light quarks completely in a first ap- 
proach and to focus on the pure glue part of the brown 
muck. However, the formalism described above is so gen- 
eral that one could also use it in a LF calculation that 
includes (dynamical) light quarks. 
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